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Abstract 

The paper is devoted to the calculation of additional two-loop 0{a 2 s ) MSSM cor- 
rections to the relation between the pole mass of the t-quark and its running mass 
in the DR scheme. Firstly, a contribution from axial part of the quark self-energy, 
which was erroneously omitted in our earlier work, is obtained. Secondly, the value 
of the second order contribution from large mass expansion in mt/MsusY is stud- 
j^. ■ ied. Finally, two-loop anomalous dimension of the running quark DR-mass in the 

O supersymmetric QCD is calculated. 
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1 Introduction 

M. 

Two-loop 0(a 2 ) MSSM corrections to the relation between the pole mass of 
the t-quark and its running DR mass had been calculated in [1] by means of 
the large mass expansion [2] in small parameter m t /Mh ar( i, where Mh ar< i stands 
for all mass scales involved in the problem that are much larger than m t . 

The initial idea of this paper was to provide more terms in this expansion and 
to study the influence of these terms on the final result. We restricted ourselves 
to the terms 0{m 2 / M^ ard ). Contrary to our initial expectations it was found 
that these terms are negligible (they affect the result of [1] approximately by 
0.1%). 
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In our previous calculations [3] we have reproduced the result obtained in [1] 
and it was found that the non-zero axial contribution to quark self-energy had 
not been taken into account by the authors of [1]. This inconsistency has been 
fixed in this paper, but as numerical analysis showed axial contribution is also 
negligible. Nevertheless, the result for the relation between the pole mass of 
the t-quark and its running mass in the DR scheme presented here is free of 
these errors. 

Another issue studied in this paper is the calculation of two-loop anomalous 
dimension of the running quark mass in the supersymmetric QCD (a subset of 
MSSM, relevant to the calculation of a 2 corrections). It was found that con- 
trary to the non-supersymmetric case the bare mass of a quark considered in 
the QCD sector of the MSSM can be used for extraction of anomalous dimen- 
sion of the quark mass, i.e. assuming its independence on renormalization scale 
p, and differentiating it with respect to log/2, we acquire correct expression for 
the two-loop anomalous dimension. 



2 The pole mass of the t-quark 

The pole mass of a particle is defined as a real part of the complex pole of the 
resumed propagator (we discuss only perturbative effects). The full connected 
propagator of a quark can be written as 



p-m- E(p, mj' 

where 

E(p, rrii) = pE y (p 2 , m 2 ) + p7 5 E A (p 2 , m 2 ) + mE s (p 2 , m,) (2) 

is the self-energy of the quark, so the pole mass M po i e satisfies the following 
equation: 

(l + E y (M^ e , m 2 )) 2 M 2 pole - E 2 (M 2 pole , m 2 )M 2 pole 

-m 2 (l-E s (M 2 oie ,m 4 )) 2 = 0. (3) 

Solving this equation perturbatively, one gets 
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M P° le m = aM' 1 » + a¥ 2 ', where (4) 
m 

M (1) = Z { v\m 2 ,m1) + Z ( s\m 2 ,m i ), (5) 
= Y%\m 2 ,mt) + ^\m 2 , mi ) + ^f(m 2 ,m 2 ) 

+M« (e£V,™?) + 2m 2 A ( S «(p 2 ,m 2 ) + 4V,™*)) l P w) (6) 

and a stands for all couplings of the theory. Using eq. (4), one calculates the 
relation between pole and running masses of the t-quark. In our case mass 
parameter m corresponds to the running mass m t (p) defined in the modified 
DR scheme [4]. In the MSSM strong interactions distinguish left- and right- 
handed particles so at the two-loop order we have to take into account the 
axial part of quark self-energy E^. A quantity we want to compute is defined 
in the following way 

Awtf = Mf ole - mt(p) 

The lagrangian of the supersymmetric QCD and relevant diagrams that con- 
tribute to the quark self-energy can be found in [1]. To evaluate these diagrams 
the large mass expansion has been used. According to its prescription asymp- 
totic expansion of a Feynman integral F r which depends on the large masses 
Mi, M 2 , . . . , small masses m 1? m 2 , . . . , and small external momenta pi,P2, ■ ■ ■ 
can be expressed as follows [2]: 

F r (pi,...,Mi,...,mi,...) = J2 F r/ 1 (p,m)M, / (p^,m)F 1 (M,m,p, / ), (8) 



where the operator A4 7 (p 7 ,m) performs Taylor expansion in small external 
(with respect to the subgraph 7) momenta and masses. The sum runs over 
all asymptotically irreducible subgraphs of the original graph T. The reason 
of such a complication is that naive expansion of a Feynman integral in small 
parameters produces spurious IR divergences which have to be subtracted by 
adding proper counter-term diagrams. 



Using two simple facts 



(1) There can be only even number of superparticle lines in a single vertex, 
since MSSM lagrangian is R-invariant, 

(2) All superparticles are considered as heavy in our problem, 

one can prove that only three types of subgraphs are possible in our problem: 

(1) All lines of a diagram are hard (Fig. 1). 

(2) A diagram with one light line. All other lines are hard (Fig. 2). 
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r 7 = r 



Fig. 1. The diagram is composed of hard lines only; therefore, the naive expansion 
can be used. 

(3) A diagram has one cut composed of two light lines 1 . All other lines are 
hard (Fig. 3). 

Expression (7) can be written as a Taylor series in mt/Mhard'- 



oo oo 



Am t , v , v -- ( „) 



m t 



i=\ n=—l 



oo oo 

= 1 + « s E <^ + «2 E m >2 ] + 0{a% (9) 

n=— 1 n=—l 

where M hard is absorbed into . We restrict ourselves to 0{m^/M^ ard ) terms, 
thus, 



1 ^ (n) , 2 M /m\ 

« 1 + a s E 0" ; + a ^ E °"2 10 



where cij 



(n) _ _ n -(n) 



We performed the calculation of (10) in a semi-automatic fashion. First, Fey- 
nArts [6] is used to generate the diagrams. Then the large mass expansion of 
individual diagrams is done by means of the C++ library prop2exp [7], based 
on GiNaC [8]. The prop2exp library performs large mass expansion according 
to (8), thus, calculation is reduced to evaluation of 2- loop vacuum integrals 
and 1-loop on-shell propagator type integrals. This task is done by the bub- 
blesll [9] C++ library, which recursively reduces 2-loop vacuum integrals to a 
master- integral [10] by integration by parts method [11]. All renormalization 
constants needed for acquiring finite answer can be found in [1]. 



3 The running mass of the t-quark and its anomalous dimension 



As was mentioned above running DR-mass of a quark is defined as renormal- 
ized quark mass in modified DR renormalization scheme. In supersymmetric 

1 expansion of such a diagram was discussed in detail in [5] 
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Fig. 2. The diagram has one light line, as a consequence, non-trivial subgraph 71 is 
also needed. 

1 







Fig. 3. The diagram has one cut composed of two light lines, so non-trivial subgraph 
71 has to be taken into account. 

QCD we have the following relation between bare quark mass m and its 
running DR-mass m(p) at the scale Jx: 



m = m(fx) {l + 



a 



2 f8Z^ sz^y 



47T £ 



(47T) 



where a s = ^/(47r), 



;n) 



5Z£ 

SZ (2,2) 



-2C F , 

-QC F + 3C A C F + 2C 2 F 



(12) 
(13) 
(14) 



and Cf = 4/3, Ca = 3 are casimirs of SU(3). Anomalous dimension of the 
quark mass 7 m is denned as 



d 



d log fi' 



;m(ii) = 7 m m(/x) 



(15) 



We consider two ways of obtaining quark mass anomalous dimension. As we 
know physical quantities do not depend on the scale parameter p, so it is 
possible to find perturbative expansion of 7 m in coupling constants of the 
theory by differentiation of pole mass expressed in terms of running parameters 
(4). Also it is possible to extract mass anomalous dimension from bare mass 
assuming that it is renormalization group-invariant. In [4] it was noticed that 
in non-supersymmetric QCD in modified DR scheme suggested by the authors 
of that work, these two definitions do not produce the same result. In this 
paper we checked that in supersymmetric QCD two procedures mentioned 
above renders the same result in spite of the fact that we had used the same 
renormalization prescription as in [4]. The reason of such a coincidence turns 
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to be in the fact that in supersymmetric QCD gauge coupling and so-called 
evanescent coupling of e-scalars renormalize in the same way, so we can set 
them equal to each other at any scale Jjl. 



The final result for anomalous dimension of the quark mass in supersymmetric 
QCD in DR looks like 



In a general theory with spontaneous gauge symmetry breaking Yukawa beta- 
functions do not coincide with anomalous dimension of corresponding fermion 
mass, see e.g. [12]. However, since the tree-level MSSM Higgs potential does 
not depend on strong coupling constant a s , and no loop corrections give con- 
tributions proportional to a s only (intuitively: the vacuum is colorless), the 
a™ coefficients in the perturbative expansion of (16) should be identical to the 
corresponding terms of top quark Yukawa coupling beta-function. Comparing 
(16) with two-loop MSSM top quark Yukawa coupling beta-function [13,14], 
one can see that a 2 s term of that beta-function indeed coincides with the cor- 
responding term of the anomalous dimension (16). This gives an additional 
confirmation of the correctness of our result. 



4 Numerical results 

The analytical result of our calculation is complicated due to the presence of a 
large number of masses, and no phenomenologically acceptable limit seems to 
exist. Therefore, we present here the numerical analysis of our results 2 . First, 
let us discuss the relative value of second order terms of large mass expansion, 
i.e. 



We investigated this quantity in the following regions of the CMSSM param- 
eter space 




(16) 



rel 



a^\M z ) + af\M z ) 



(17) 




600 GeV <m < 1000 GeV 
200 GeV < mi/2 < 1000 GeV 



2 The full answer is available at http://theor.jinr.ru/~varg/dist 



in a form of C++ 



library 
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Fig. 4. Relative value of second-order term of large mass expansion a 2 el (see (17)) 
as a function of m-i/2- 

for ao = 0, 0.4mo, — 0.4m and tan/5 = 33, 50. We only consider fi > and 
large values of tan (3, since small tan (5 and negative /x seem to be excluded by 
experimental data [15]. We found that a 2 cl < 10~ 3 in these regions. Typical 
behaviour of a 2 cl is shown in Fig. 4. Thus, contribution of second order large 
mass expansion terms to the relation between the pole and running masses 
(7) is negligible. 

Numerical values of running SUSY parameters at the M z scale have been cal- 
culated as a function of CMSSM parameters with heavily modified version of 
SOFTSUSY [16] in the framework of the mSUGRA supersymmetry breaking 
scenario. 



In the original SOFTSUSY code (as of version 1.9) two-loop SQCD contri- 
bution to the relation (7) is neglected, while two-loop QCD contribution to 
that relation is taken into account. This approximation is not applicable in 
the region of CMSSM parameters we considered. Figure 5 demonstrates this. 
Fig. 5 shows my 2- an d m -dependence of the 0{a 2 s ) supersymmetric QCD 
corrections to (7) (here and in what follows all SQCD corrections include con- 
tributions from "pure QCD" diagrams). For comparison we plotted 1-loop 
SQCD contribution and 2-loop QCD contribution. One can see that 

(1) two-loop SQCD correction is about of 30% of the one-loop one, 
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Fig. 5. Two-loop (af) and one-loop (a s ) SQCD contributions to 
Amt{Mz)/mt(Mz)- Solid lines correspond to SQCD corrections, dashed line 
correspond to "pupe QCD" ones. 
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(2) two-loop SQCD contribution is of the same order of magnitude as one- 
loop QCD correction, 

(3) two-loop QCD contribution is about of 20% of total two-loop SQCD 
correction. 

Thus, two-loop SQCD correction to the relation between the pole and running 
masses of the t-quark is not negligible and should be taken into account in 
phenomenological analysis of MSSM 3 . 

Taking into account this correction yields sizable change (more than 10%) 
of predicted masses of heavy Higgs bosons and chargino (see Fig. 6). This 
fact represents strong dependence of \x and Ma on the heavy quark Yukawa 
couplings in certain regions of parameter space [18]. This dependence leads to 
relatively large discrepancies [18,19] between the values of predicted masses 
given by different software for calculation of MSSM mass spectrum [16,17]. It 
should be noted that change due to two-loop SQCD correction to the relation 
(7) exceeds these discrepancies. 

However, masses of squarks, gluino, and relatively light particles (lightest neu- 
tralino, lightest Higgs boson) do not obtain any significant changes due to 
above mentioned two-loop SQCD contribution, see Fig. 7. 



5 Conclusion 



In this paper, we presented the result of calculation of the two-loop corrections 
to the relation between pole and running masses of the t-quark in the super- 
symmetric QCD. We provided a numerical analysis of the value of these correc- 
tions in different regions of the CMSSM parameter space. Our analysis showed 
that calculated second-order terms of large mass expansion 0(m^/M^ ard ) give 
negligible contribution to this relation, so zero-order result given in [1] provides 
reliable approximation for phenomenological studies of the MSSM. Analysis 
given in [18] demonstrates that there exists two regions of the MSSM param- 
eter space where accurate predictions based on computer codes [16,17] are 
difficult: the large tan/3 and focus-point regimes. These two regions require a 
more precise determination of the heavy quark Yukawa couplings, or equiva- 
lently a more precise determination of running quark masses. We showed that 
two-loop supersymmetric QCD corrections give sizable contribution to (10) 
and has to be included in computer codes used to calculate MSSM spectra 
[16,17]. 



3 Alternatevly, if one neglects 2-loop SQCD correction, one can neglect 2-loop QCD 
correction as well 
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200 400 600 800 1000 200 400 600 800 1000 

mi/2, GeV m 1/2 , GeV 



Fig. 6. Masses of heavy Higgs bosons and chargino as a functions of mj^. Dashed 
lines correspond to the output of the original SOFTSUSY code (two-loop SQCD 
contribution to Am t (Mz) jmtiMz) is neglected), solid lines correspond to the same 
code with above mentioned SQCD corrections. Gray regions represent discrepancies 
between the values of masses predicted by different programs for calculation of 
MSSM mass spectrum [19]. Here tan/3 = 50, A = 0, m = 1000 GeV 



As a by-product of our calculation, we also obtained two-loop anomalous di- 
mension of the running quark DR-mass in the supersymmetric QCD. 
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Fig. 7. Change of predicted mass of the lightest Higgs boson and the lightest neu- 
tralino mass due to two-loop SQCD contribution to Amt(Mz) /m t (Mz) as a func- 
tion of mi/2- Here mo = 1000 GeV, tan/3 = 50, Aq = 0. 
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